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ABSTRACT 

This  paper  deals  with  the  blind  identification  of  the  vari¬ 
ances  of  the  additive  noises  in  a  two-channel  SIMO  system. 
They  are  estimated  by  using  a  frequency  domain  errors-in- 
variables  scheme,  when  finite  impulse  responses  (FIR)  filters 
are  driven  by  a  white  input  sequence.  It  should  be  noted 
that  the  proposed  approach  has  the  advantage  of  addressing 
the  unbalanced  case,  i.e.  when  the  noise  variances  are  not 
the  same.  Moreover,  the  algorithm  enables  the  solution  to 
be  directly  obtained  from  the  noisy  data  and  no  optimization 
algorithm  is  required. 

1.  INTRODUCTION 

For  the  last  years,  a  great  deal  of  interest  has  been  paid  to  sin¬ 
gle  input  multiple  output  (SIMO)  system  identification  [1]- 
[5]  in  signal  processing.  Thus,  in  the  framework  of  speech 
enhancement  based  on  two  microphones,  blind  identifica¬ 
tion  methods  such  as  [6]  have  been  proposed  to  retrieve  the 
speech  signal  from  reverberated  observations  disturbed  by 
additive  noises.  In  mobile  communications,  identification 
approaches  have  also  been  of  interest  to  carry  out  equaliza¬ 
tion  of  the  received  signal  [7]. 

In  the  following,  let  us  consider  the  SIMO  system  de¬ 
picted  in  figure  1,  where  s(k)  denotes  the  input  signal  of 
M  =  2  filters,  defined  by  their  finite  impulse  responses  (FIRs) 
hi(k),  ( i  =  1,2),  with  maximal  order  L.  The  filter  outputs 
Xj(k)  are  then  disturbed  by  additive  zero-mean  noises  bfk), 
with  variance  of,  (i  =  1,2).  This  yields  the  noisy  observa¬ 
tions  yfk),  ( i  =  1,2).  Moreover,  s(k)  is  assumed  to  be  a  sta¬ 
tionary  zero-mean  white  process  with  unit-variance,  without 
loss  of  generality. 


cross-relations  between  the  filtered  outputs: 

hj(k)  *x,(k)  =  hj(k)  *  (hfik)  *  s(k)) 

=  ( hi(k)*hj(k ))  *s(k)  =  hj(k )  *Xj(k).  (1) 

For  this  purpose,  they  search  the  null  space  of  a  suit¬ 
able  positive  hermitian  block  Toeplitz  matrix,  built 
with  the  cross-correlation  coefficients  of  the  noisy  fil¬ 
ter  outputs.  As  an  alternative,  the  subspace  based 
method  presented  in  [2]  consists  in  retrieving  the 
span  of  the  autocorrelation  matrix  R/  of  the  vector 
x(&)  =  [x[ (k)  ■  ■■xrM(k)\ T  with  xfik)  =  [xfik)  ■  ■  ■xi(k-L)]r . 
It  should  be  noted  that  both  methods  are  equivalent  when 
dealing  with  two  channels. 

Adaptive  approaches  have  also  been  studied,  either  in 
the  time  domain  [3]  or  in  the  frequency  domain  [4]  [5], 
They  aim  at  minimizing  the  cross-relations  error,  i.e. 
E\\hj(k)  *y,(k)  — hfk )  *yj(k) |2]  and  are  based  on  the  least 
mean  square  (LMS)  and  Newton  algorithms  applied  to  the 
noisy  data. 

However,  the  above  methods  cannot  be  used  when  the  ad¬ 
ditive  noises  do  not  have  the  same  variance,  i.e.  when  deal¬ 
ing  with  the  “unbalanced"  case.  Few  methods  deal  with  this 
issue.  In  [6]  and  [8],  the  approach  is  based  on  the  Frisch 
scheme  [9],  which  is  a  class  of  errors-in-variables  (EIV)  so¬ 
lutions.  The  EIV  models  [11],  initially  developed  in  the  fields 
of  statistics  and  identification,  assume  that  the  available  data 
are  disturbed  by  additive  error  terms.  Given  a  generic  process 
described  by  K  variables  {tk}k=i,-  ,k,  the  formulation  of  an 
EIV  estimation  problem  consists  in  determining,  on  the  only 
basis  of  noisy  observations  {za  =  4 +  «a}a=i, •••.*:-  the  set  of 
X-tuple  {fk}k=i,-  -,K  that  satisfies: 


bx(k),  (of) 


Figure  1:  two-channel  SIMO  system 

SIMO  system  identification  consists  in  estimating  the 
FIRs  and  may  include  the  estimations  of  the  variances  of  the 
additive  noises1.  Thus,  Xu  et  al.  [1]  exploit  the  so-called 

1  It  should  be  noted  that,  in  some  specific  cases,  the  variance  of  each 


X\t\  +  X2h  5 - b  Xmk  =  0 

or  equivalently, 

R[  ['ll  X2  ■■■  Xk]T  =0  (2) 

where  Rt  is  the  covariance  matrix  of  {tk}k=i,  -,K  and  the  up- 
perscript 7  denotes  the  transpose. 

When  each  noise  term  tik  is  assumed  to  be  zero-mean,  inde¬ 
pendent  of  every  other  noise  term  and  every  variable  4,  one 
has: 

Rt  -  Rz  Rn 

where  Rz  and  Rn  respectively  denote  the  covariance  matrix 
of  {zk}k=\,  - ,K  and  of  {ha}a=i, 

additive  noise  could  be  estimated  independently  of  each  other.  In  speech 
processing,  it  could  be  updated  during  silent  frames.  Subspace  methods 
[10]  could  also  be  considered  provided  that  the  spectrum  of  the  noise-free 
signal  is  discrete.  However,  these  cases  are  too  restrictive  and  we  propose 
here  to  study  a  more  general  solution. 
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At  that  stage,  the  so-called  Frisch  scheme  [9]  consists 
in  searching  for  the  diagonal  matrices  Rn  which  enable 
Rz  —  Rn  to  be  positive  semidefinite: 

Rz  Rn  >  0  (3) 

When  dealing  with  SIMO  system  identification,  given  the 
cross  relations  (1)  and  by  selecting  K  =  2L  +  2,  the  set  of 
2L  +  2  variables  {tk}k=i,..,2L+2  can  be  for  instance  defined  as 
follows: 

[t\---tiL+2]T  =  [x2(k)---x2{k-L)  -XI {k) - XI (k-L)]T , 

whereas  the  (2L  +  2)-tuple  (Ai ,  -  -  -  ,^2^+2)  satisfies: 

[A,  ■  •  -Wf  =  [hi{0)  ■  ••Ai(L)  h2{ 0)  •  •  -h2{L)]T . 


Taking  into  account  (4)  and  (5),  the  noise  variances  are  solu¬ 
tions  of  the  following  equation  system: 


X<xua2{w)  :=  det 


cq  0 
0  a2 


=  0 


trace 


cq  0 
0  a2 


>0  (w€  [ — 7r; 7r]) 


(6) 

(7) 


where  oq  and  a2  are  the  unknowns. 

Estimating  the  noise  variances  is  very  useful  as  it  makes 
it  possible  to  estimate  X  (w)  and  consequently  to  retrieve 
the  impulse  responses  hj(k')  ( i  =1,2)  since  equation  (2),  or 
equivalently  equation  (1),  can  be  rewritten  in  the  frequency 
domain  as  follows: 


In  that  case,  Zk  corresponds  to  the  noisy  observations.  Ac¬ 
cording  to  [8],  solving  (3)  yields  an  infinite  set  Si  of  vari¬ 
ances  (csf, of).  In  [6],  we  suggest  carrying  out  once  again 
the  Frisch  scheme  by  selecting  K  =  21  +  2  with  l  >  L  to  find 
the  true  variances  which  belong  to  Sl  D  S/.  Then,  the  impulse 
responses  can  be  retrieved  according  to  (2).  However,  there 
are  two  main  drawbacks: 

1.  There  may  be  spurious  solutions  if  /  is  not  suitably  cho¬ 
sen  ( l  —L  must  not  be  too  low). 

2.  We  do  not  take  benefit  of  the  fact  that  the  true  variances 
belong  to  every  set  Si  with  l  >  L.  Indeed,  for  the  sake  of 
simplicity  and  computation,  only  two  Frisch  schemes  are 
carried  out. 

Therefore,  in  this  paper,  we  propose  to  develop  a  frequency 
domain  EIV  approach  for  SIMO  identification.  This  method 
avoids  the  above  drawbacks.  The  rest  of  the  paper  is  orga¬ 
nized  as  follows:  section  2.1  deals  with  the  variance  estima¬ 
tion  issue.  The  identifiability  conditions  and  the  proposed 
algorithms  are  introduced  in  section  2.3.  Then,  experimental 
results  are  presented  in  section  3. 


Ker(AH(w))  =  Span 


H 1  (w) 

H2(w) 


In  the  next  two  subsections,  we  will  study  how  to  solve 
the  system  (6)-(7),  and  will  present  algorithms  for  the  noise 
variance  estimation. 


2.2  First  approach  to  blindly  estimate  the  variances 

In  addition  to  the  trivial  solution  (of.  of),  the  system  (6)  and 
(7)  admits  another  solution  if  and  only  if  one  can  find 

(y,p)  €  Rf  x  R*  such  as 

Xn{w)-yX22(w)=  p.  (8) 


See  appendix  1  for  proof. 

Then,  a  first  approach  we  could  consider  to  blindly  es¬ 
timate  the  noise  variances  would  consist  in  minimizing  the 
following  cost  function: 


Ji(aha2) 


X, 


Ct\  ,0^2 


(w) 


=  (X, 


,a2  v 


X, 


a[,a2 


(9) 


2.  SIMO  SYSTEM  IDENTIFICATION  IN  THE 
FREQUENCY  DOMAIN 


where  <•!•>  denotes  the  hermitian  product  defined  as  fol¬ 
lows: 


2.1  Focus  on  the  variance  estimation  issue 


Given  the  assumptions  on  the  input  signal  s(k)  made  in  the 
introduction  above,  the  power  spectrum  density  matrix  of  the 
vector 

x(k)  =  [x'2 {k)  — xi  (k)]T  is  defined  by: 


X(w) 


Xu  (w)  X\2{w) 

X2\ (w)  X22(w) 

\H2(w)\2  — H2(w)H^(w ) 

_—H* (w)H2(w)  \H\  (w)|2 


X(w)  is  positive.  In  addition,  it  is  singular  since  one  has 


det(Z(w))  =  0. 

Both  properties  are  therefore  denoted  as  follows: 


X(w)>  0. 


(4) 


However,  only  the  power  spectrum  density  Y  (w)  of  the  noisy 
observation  vector  y(k)  =  \y2(k)  —  yi(&)r  is  available. 
According  to  figure  1,  it  satisfies: 


Y(w) 


Tn(w)  Yi2(w) 

l2i(w)  Y22(w) 


X(w) 


Oy  0 

0"  of 


X(w)  +B(w) 


(5) 


(A(w)\B(w))  =  —  [  A(w)B*(w)dw.  (10) 
1  2 n  J-x 

The  hermitian  product  (10)  can  be  expressed  by  using 
the  Fourier  coefficients  of  the  2^-periodic  functions  A{w) 
and  B(w),  denoted  ak  and  bk  respectively.  Indeed,  since 
A(yv)  =  ^  akejk"  and  B(w)  =  bkekkw,  one  obtains 

ke  z  Ar6Z 


<A(w)  |B(w)>  =  i-  f  (  £  ak^)  (  £  fcV'M *dw 


1 

2  % 


p7l 

/  T 

J  K  n€Z 

=  £ak(bk)*. 

ke  z 


Eak(bk~ny 

kez 


ejmdw 


Therefore,  when  dealing  with  hermitian  products  between 
the  elements  of  the  power  spectrum  density  matrices  such 
as  Xjj(w)  for  instance  (i,j  =  1,2),  the  hermitian  product  can 
be  hence  defined  from  the  intercorrelation  sequences  of  xt(k) 
and  Xj(k). 

The  gradient  of  J\  (cq ,  a2)  is  given  by 


dJ{aua2)  / dXai,cc2(w )  v  /  x\  z'-  1  ox 

Ba,  =2\  ila,  0  =  1.2). 
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When  cancelling  it,  the  following  two  equations  have  to  be 
solved: 


6 


Ji  (01,0:2)  theorique 


«1  =  02  («2)  (11) 

y°i  i  af  -  (g°-f-  (F22  (w)  |  Fi  i  (w) ))  a2  +  (*22  (w)  |  G(  w) ) 
a2  -2y^2a2  +  ||F22(w)||2 

«2  =  01  («l)  (12) 

^“t  -  (^°+ (y!  1 W  1*22 (w) ))  oq  +  (Fi  1  (w) I G(w) ) 

a\  ~2y|1ai  +  ||Fn(w)||2 

The  solution  of(ll)end(12)  are  given  by  searching  the  roots 
of  a  5th  degree  real  polynomial,  depending  on  either  oq  or  a2. 


In  theory,  the  variances  of  the  additive  noises  can  be  ob¬ 
tained  as  follows: 

•  when  the  trivial  solution  (of,  of)  is  the  only  solution  to 
the  equation  (6),  it  can  be  obtained  by  minimizing  the 
cost  function  J\  (oq ,  a2). 

•  When  there  are  two  solutions  to  the  equation  (6),  that 
can  be  also  found  by  minimizing  J\  (oq ,  a2),  there  exists 
y  e  R*  and  peR*  such  as 

Xn{w)  -  yX22(w)  =  p .  (13) 

See  appendix  2.  However,  when  7  >  0,  only  the  trivial 
solution  satisfies  (7). 

Nevertheless,  in  real  cases  and  when  (6)  admits  two  so¬ 
lutions,  the  cost  function  J\(tt\,  0:2)  can  then  have  very  small 
values  in  the  neighborhood  of  theses  solutions.  As  a  conse¬ 
quence,  retrieving  them  becomes  a  hard  task,  due  to  the  er¬ 
rors  introduced  when  computing  Y  (w)  from  a  finite  number 
of  available  samples: 

•  on  the  one  hand,  B(w)  is  not  diagonal. 

•  On  the  other  hand,  the  unbiased  estimates  of  the  intercor¬ 
relation  sequences  are  considered. 

To  illustrate  this  point,  let  us  consider  the  example  given  be¬ 
low  where  the  two  following  FIRs  satisfy: 

M  =  {2, 1,-1, 2} 

,  (4  -65  +  V557A  65-V552A  t 

*2  =  1 3  ’ - 54 - ' - 24 - 'T 


In  figure  2  the  theoritical  cost  function  /i(ai,a2)  plot¬ 
ted  in  dB  exhibits  the  trivial  solution  and  the  spurious  one. 
Nevertheless,  it  takes  very  small  values  in  the  neighborhood 
of  both  solutions.  In  figure  3,  when  a  real  simulation  is  com¬ 
pleted  with  N  =  10000  observations  samples  of  y  1  and  y2,  the 
estimated  cost  function  /1  (ai ,  a2)  does  not  make  it  possible 
to  point  out  the  true  solution. 

More  generally,  we  will  see  that  this  problem  happens  when 


(Xn(w)|X22(w))| 


Xn(w) 


X2l(w) 


where  An  (w)  abd  A22(w)  are  defined  as  follows: 

Aii(w)  =An(w)+Xii  (14) 

A22(w)  =A22(w)+x?2,  (15) 

with  Jtjj  =<  A{j(w)|l  >  and  x\2  =<  1 1  >  denoting 

the  mean  values  of  Ai  1  (w)  and  Aj  1  (w)  respectively.  This  will 
be  confirmed  in  the  next  subsection. 


Figure  2:  theoritical  criterion  J\  (oq ,  a2)  when  a  spurious  so¬ 
lution  to  (6)  exists.  Here,  7  >  0. 

J\  (0:1,02),  a  real  case 


Figure  3:  criterion  /i(ai,a2)  in  dB,  in  a  real  case  while  a 
spurious  solution  to  (6)  should  have  existed.  N  =  10000. 


2.3  Second  approach  to  blindly  estimate  the  noise  vari¬ 
ances 

We  propose  a  new  criterion  to  estimate  the  noise  variances 
when  the  equation  (6)  has  a  single  solution.  It  is  based  on  the 
following  property.  Let  us  assume  that  both  Ai  1  (w)  /  0  and 
2f22(w)  ^  0,  and  An(vv)  ^  A22.  According  to  appendix  2,  we 
have: 


(a)  XauOl2(w)=0  has  a  single  solution 

( b )  Acq  (w)  =  0  has  a  single  solution  <£=> 


(c)  |(Ah(w)|A22(w))|  <  An  (w)  X22  (w) 


(16) 


where  A ai,a2(w)  is  defined  as  in  (14)-(15). 

Then,  when  there  is  a  single  solution  to  equation  (6),  it  can 
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be  found  by  solving  equation 


Algorithm  1 


Algorithm  2 


-^0:1,02  (w)  — 

Therefore,  we  propose  to  minimize  the  following  criterion: 


J2(aua2)  = 


X, 


a  [,<X2 


O) 


G(w)  -  a\%2  (w)  -  a2Ti  i  (w)+a\  a2 


(17) 


where  G(w)  =  Fu(w)F22(w)  —  F2i  (w)Fi2(w)  and  G(w)  is  de¬ 
fined  as  in  ( 14)-(  1 5). 

J2(a\,a2)  is  a  quadratic  function.  Equating  to  zero  the  gra¬ 
dient  of  J2(a i ,  a2)  leads  to  the  following  equations: 


—  F22 (w)  G(w)  -  a2Tn(w)  -  aiF22(w)f  =  0 
-Fii(w)  G(w)  -  ai?22(w)  -  a2?n(w)\  =  0 


One  finally  deduces: 

-  (?22  (w)  I  ?! !  (w) )  a2  ±  (?22  (w)  I  G(w) ) 

oq  =  — - - =4 - - - -  (18) 

II^22(W)||^ 

-<T22(w)|Tii(w))a1  +  (Fi1(w)|G(w)) 

C&2  =  - ~ -  (IV) 

IMI2 

Hence,  this  linear  system  can  be  re-written  as  follows: 


(?22(M')|?ll(w)) 

ifeM!!2 

(y22(w)\G(w)) 

0=1 

P22  W!|2 

(  y22  (w)  |  iq  i  (w) ) 

IIHiMlP 

a2 

(rn(w>)|G(w)) 

-  HEiMII2  . 

■v- 

M 


One  retrieves  the  condition  (16)-(c)  by  imposing  det(M)  >  0. 

It  should  be  noted  that  the  condition  (7)  is  no  longer  nec¬ 
essary  when  equation  (6)  is  assumed  to  have  a  unique  solu¬ 
tion. 


3.  RESULTS 

We  propose  to  illustrate  the  algorithm  with  two  examples. 

The  first  example  deals  with  a  L  =  2nd  order  system,  the 
FIRs  corresponding  to  the  following  transfert  functions  are 
defined  by: 

H\(z)  =  1  ±z±z_1  and  H2(z)  =  1  -z  +  z~1- 

N  =  150  samples  are  assumed  to  be  available.  Moreover,  the 
two  curves  in  the  plan  (oq ,  a2)  defined  by  the  relations  (11) 
and  (12)  for  the  first  algorithm,  and  the  two  lines  in  the  plane 
(o=i,  a2)  defined  by  the  relations  (18)  and  (19)  for  the  second 
algorithm,  are  represented  in  figure  4. 

In  the  second  example,  the  proposed  algorithms  are 
then  tested  with  higher-order  systems  (L  =  100),  whose  time 
and  frequential  domain  representations  are  given  in  figure  5. 
Here,  N  =  1000  samples  are  available.  The  good  perfor¬ 
mances  are  illustrated  in  figure  6. 

Numerical  results  based  on  Monte  Carlo  simulations  are 
given  in  Table  1  and  2,  respectively  for  the  examples  1  and  2. 
In  each  case,  the  signal-to-noise  ratios  (SNR)  for  each  chan¬ 
nel  are  respectively  equal  to  SNR \  =  10 dB  and  SNR2=5dB. 
The  estimations  of  the  variances  provided  by  the  cost  func¬ 
tions  7i(ai,a2)  (9)  and  /2(ai,a2)  (17)  are  compared.  In 
every  case,  only  few  samples  (i.e.  a  few  hundred)  are  re¬ 
quired  to  estimate  the  pair  (of,  )• 


Figure  4:  noise  variances  estimation  for  example  1. 


hi(k)  h2(k ) 


Figure  5:  L,h  order  FIR  and  their  spectral  representations. 


Algorithm  1 


Algorithm  2 


Figure  6:  noise  variances  estimation  for  example  2. 


Table  1:  variance  estimation  and  standard  deviation  based  on 
100  noise  realizations,  for  the  SIMO  system  of  example  1. 


example  1 

erf  =  0.290 

of  =  0.998 

Criterion  J\  [a\ ,  a2) 

0.383  ±0.054 

0.930  ±0.100 

Criterion  J2{a\ ,  a2) 

0.498  ±0.121 

1.042  ±0.173 

Table  2:  variance  estimation  and  standard  deviation  based  on 
100  noise  realizations,  for  the  SIMO  system  of  example  2. 


example  2 

erf  =  0.776 

of  =2.318 

Criterion  J\  ( oq ,  a2) 

1.252  ±0.357 

1 .836  ±  0.338 

Criterion  J2{a  \ ,  a2) 

1.506  ±0.447 

2.188  ±0.427 

4.  CONCLUSION 

In  this  paper,  the  variance  estimation  issue  in  the  unbalanced 
case  is  addressed.  For  this  purpose,  we  propose  a  frequency 
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domain  errors-in-variables  approach.  More  particularly,  we 
investigate  a  first  criterion  and  its  the  limits  and  hence  pro¬ 
pose  an  alternative  one.  This  criterion  has  the  advantage  of 
being  simpler  than  the  previous  one  since  it  is  quadratic. 

5.  APPENDIX  1 

Let  us  assume  that  ((*1,052)  7^  (of,C rf)  is  a  solution  of  (6). 
Then,  ai  f  cry  AND  0.7  f  of.  By  rewritting  (6)  as  follows: 

(Xn(w)  +  al  -  a^(x22(w)  +  af  -  (*2)  =  XU(w)X22(w), 

there  is  necessarily  y  £  R*  such  that: 

Xn(w)  +  -  ai  =  yX22 (w)  (20) 

,  1 

X22(w)  +  ol  —  a2  = -Xu(w).  (21) 

Moreover,  since  the  pair  (051,0:2)  also  satisfies  (7),  one  de¬ 
duces: 


since  both  Xn(vv)  7^  0  and  X22(w)  f=0.  As  a  consequence, 
there  is  y  S  W  such  as: 

(of  -  ai)  =  -y(of  -  a2)  (23) 

Xu  =  yX22.  (24) 

Substituing  (23)  and  (24)  in  (22 )-(ii)  yields  the  following  re¬ 
lation: 

ai  =  of  -  yt22  +4t  (25) 

052  =  of -^-+4,  (26) 

Finally,  since  Xn(vv)  f  X22,  —  yrfj  +4i  7^  0-  the  pair 
(cq ,  052 )  defined  by  the  relations  (25)  and  (26)  satisfies 
G!i  /  of  and  0.7  /  of,  and  is  solution  of  (6),  that  contradicts 
(a). 

(b)  •<=>■  ( c )  since  both  (£>)  and  (c)  means  that  Xi  1  and  X22 
are  not  colinear. 


X\l(w)  +  X22  (vv)  +  of  -  051  +  of  -  052 
=yX22(w)  +  ixn  (w)  >  0 

Necessarily,  y  >  0.  By  introducing  p  =  ai  —  O?  and  in  the 
relation  (20),  one  obtains  (8). 

Now,  one  assumes  thatXn(w)  —  yX72(w)  =  p  with 
y  >  0  and  pgR*.  Then,  the  pair  (ai,a2)  =  p  +  of ,  of  —  £ 
is  a  non-trivial  solution  of  (6).  Moreover,  one  can  write: 


trace 


-P  0 


(Xn(w)-p)  +  (X22(H0  +  p 


According  to  relation  (8)  and  since  both  y  >  0  and  the  func¬ 
tions  X\\(w)  and  X22(w)  are  positive,  one  has: 

Xn(w)  —  p  =  yX22  (w)  >  0 

X22(w)  +  —  =  -Xn(w)  >  0 

y  y 

Finally,  one  obtains 


trace 


So,  (p  +  of ,  of  —  y )  satisfies  the  relation 


>0. 

(7). 


6.  APPENDIX  2 

Proof :  Let  us  first  expand  the  equation  (6)  as  follows: 

Xai,a2  ('4  =  0  44- 

f(i)  Xaua2(w)  =  0 

l  (if)  (of  -  052)4i  +  (of  -  ai )x°22  (22) 

[  +(of-a2)(of-a1)  =  0. 

where  Xai  (w)  can  be  developped  as  follows: 

Xaua2(w)  =  (of-a2)Xn(w)  +  (a^-ai)X22(w). 

Hence,  ( b )  =2*  (a).  To  prove  the  converse  part,  let  us  assume 
a  non-zero  solution  to  (22 )-(/).  Then  o.\  /  of  AND  ai  /  of 
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